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Abstract 

Stochastic orders are binary relations defined on probability distributions which capture intuitive 
notions like being larger or being more variable. This paper introduces stochastic ordering of interference 
distributions in large-scale networks modeled as point processes. Interference is a major performance- 
limiting factor in most wireless networks, thus it is important to characterize its statistics. Since closed- 
form results for the distribution of interference for such networks are only available in limited cases, it 
is of interest to compare network interference using stochastic orders, for two different point processes 
with different fading or path-loss scenarios between the interferers and the receiver. In this paper, 
conditions on the fading distribution and path-loss model are given to establish stochastic ordering 
between interferences. Moreover, Laplace functional ordering is defined between point processes and 
applied for comparing interference. Monte-Carlo simulations are used to supplement our analytical 
results. 
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I. Introduction 

Since interference is the main performance-limiting factor in most wireless networks, it is 
crucial to characterize its statistics. The interference mainly depends on the fading channel 
(interfering power distribution), the path-loss model (signal attenuation with distance), and 
network geometry (spatial distribution of concurrently transmitting nodes). The spatial location 
of the interferers can be modeled either deterministically or stochastically. Deterministic models 
include square, triangular, and hexagonal lattices [1], [2], which are applicable when the location 
of the nodes in the network is constrained to a regular structure. On the other hand, only a 
statistical description of the location of the nodes is available in some scenarios. In both cases, 
the locations of transmitting nodes in the network are seen as the realizations of some point 
processes [3]-[5]. For certain classes of node distributions, most notably Poisson point processes, 
and certain attenuation laws, closed-form results are available for the interference distribution 
which determine the network performance [6], [7] (and the references therein). However the 
interference distribution is not tractable in most other cases. 

Successful transmission probability in the presence of interference can be calculated by 
determining the Laplace transform of interference [5], [8], [9]. However, closed-form expressions 
for the Laplace transform of interference are not tractable in many cases. We approach this 
problem from a stochastic ordering perspective, which is a partial order on random variables 
[10], [11]. Concepts of stochastic ordering have been applied to scenarios of interest in wireless 
communications in [12]. We will use these concepts to understand the interference in large scale 
wireless networks. In [11], [13], [14], the application of this set of tools in communication 
networks can be found. In [11], [13], the stochastic ordering has been used in studying a class 
of queueing networks. Directionally convex ordering of different point processes and its integral 
shot noise fields which are inherent from the point processes has been studied in [14]. To the 
best of our knowledge, there is no study of conditions on fading channels and path-loss models 
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for stochastic ordering of network interference in the literature. In this paper, we use stochastic 
ordering theory to compare network performance under conditions of differing fading on the 
interference link, and different path-loss models for the establishment of stochastic ordering 
of interference from different point processes. Using the conditions, we compare performance 
without having to obtain closed-form results for a wide range of performance metrics. We also 
compare different point processes which are commonly used in the literature using stochastic 
orders, and advocate Laplace functional ordering of point processes over directional convex 
ordering when interferences due to these point processes are compared. 

The paper is organized as follows: In Section II, we introduce mathematical preliminaries and 
present the system model and assumptions. Section III introduces performance metrics involving 
stochastic ordering of interference. In Section IV and V, we derive the conditions on fading 
channel and path-loss model for stochastic ordering of interference respectively. Section VI 
defines Laplace functional ordering and compares different point processes. Section VII presents 
results from numerical simulations to corroborate our claims. Finally, the paper is summarized 
in Section VIII. 

II. Mathematical Preliminaries and System Model 

Here we give a brief overview of some terminology and mathematical tools and introduce the 
system model assumed in this paper. 

A. Stochastic Ordering 

First we briefly review some common stochastic orders between random variables, which can 
be found in [10], [11]. 

1 ) Usual Stochastic Ordering: Let X and Y be two random variables (RVs) such that 

P (X > x) < P(Y > x) ,-oo <x < oo. (1) 
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Then X is said to be smaller than Y in the usual stochastic order (denoted by X < st Y). 
Roughly speaking, (1) says that X is less likely than Y to take on large values. To see the 
interpretation of this in the context of wireless communications, when X and Y are distributions 
of instantaneous SNRs due to fading, (1) is a comparison of outage probabilities. Since X,Y 
are positive in this case, x G M + is sufficient in (1). 

2) Laplace Transform Ordering: Let X and Y be two non-negative random variables such 
that 

C x (s) := E[exp (sX)} > E[exp (sY)} = C Y (s) for s > 0. (2) 

Then X is said to be smaller than Y in the Laplace transform (LT) order (denoted by X < Lt Y). 
For example, when X and Y are the instantaneous SNR distributions of a fading channel, (2) can 
be interpreted as a comparison of average bit error rates for exponentially decaying instantaneous 
error rates (as in the case for differential-PSK (DPSK) modulation and Chemoff bounds for other 
modulations) [12]. The LT order X < Lt Y is equivalent to 

E[l(X)]>E[l(Y)], (3) 

for all completely monotonic (cm.) functions l(-) [11, pp. 96]. By definition, the derivatives of 
a cm. function l(x) alternate in sign: (— l) n d n l(x)/dx n > 0, for n — 0, 1, 2, . . . , and x > 0. An 
equivalent definition is that cm. functions are positive mixtures of decaying exponentials [11]. 
A similar result to (3) with a reversal in the inequality states that 

X < U Y ^>E[l(X)} <E[l(Y)}, (4) 

for all /(•) that have a completely monotonic derivative (c.m.d.). Finally, note that X < st Y =>- 
X < L t Y. This can be shown by invoking the fact that X < st Y is equivalent to E[/(X)] < 
E[/(F)] whenever /(■) is an increasing function [11], and that c.m.d. functions in (4) are increas- 
ing. 
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B. Point Processes 

Point processes have been used to model large-scale networks [3], [5], [7], [15]— [21]. In this 
paper, we focus on stationary and isotropic point processes. A point process $ is stationary if its 
distribution is invariant to translations and is isotropic if its distribution is invariant to rotations. 
In what follows, we introduce some fundamental notions that will be useful. 

1) Campbell's Theorem [3], [5]: It is often necessary to evaluate the average sum of a function 
evaluated at the point of the process $. Campbell's theorem helps in evaluating such sums. For 
any non-negative measurable function u, 



The intensity measure A of $ in (5) is a characteristic analogous to the mean of a real-valued 
random variable and defined as A(B) = E [$(B)} for bounded subsets B C R d . So A(B) is 
the mean number of points in B. If <E> is stationary then the intensity measure simplifies as 
A(£?) = X\B\ for some non-negative constant A, which is called the intensity of $, where \B\ 
denotes the d dimensional volume of B. For stationary point processes, the right side in (5) is 
equal to A J Rd u(x)dx. Therefore, any two stationary point processes with same intensity lead to 
equal average sum of a function (when the mean value exists). 

2) Laplace Functional of Point Processes [4]: The Laplace functional L of point process $ 
is defined by the following formula 



where u(-) runs over the set IX of all non-negative functions on M d . The Laplace functional 
completely characterizes the distribution of the point process. As an important example, the 
Laplace functional L of stationary Poisson point process $ PPP is 




(5) 



L&(u) = E [e-E-e*"^)] = E e -^ d " (:r) * {d:r) 



(6) 




where A is the intensity. 
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C. System Model 

As shown in Fig. 1, we assume a transmit/receive pair communicating over a wireless channel. 
The receiver is being interfered by interference sources distributed as stationary and isotropic 
point process. The point process describes all interfering nodes [18]. Both the transmitter and 
receiver are fixed and are not considered a part of the point process. The accumulated interference 
to the receiver at the origin is of interest to quantify and it is given by 



where $ denotes the set of all interfering nodes which is modeled as a point process on R d and 
hj is a positive random variable capturing the (power) fading coefficient between the receiver 
and the x th interfering node. Here, typically d = 2 or d = 3, though this assumption is not 
necessary. The spatial region containing the interferers is assumed to be an infinite area [4], [5], 
[7], [16]. Moreover, {h\ x) } x are i.i.d. random variables and independent of the point process. 
The path-loss is captured by a function g(-) : R + — > R + which is a continuous, positive, non- 
increasing function of ||x|| and assumed to depend only on the Euclidean distance ||x|| from the 
node x to the receiver at the origin. We consider the following general path-loss model [4]-[7]: 



for some b > 0, 5 > d and a e {0, 1}, where 5 is called the path-loss exponent, a determines 
whether the path-loss model belongs to a singular path-loss model (a = 0) or a non-singular 
path-loss model (a = 1), and b is a compensation parameter to keep the total receive power 
normalized which will be discussed in detail in Section V. It is noted that from Campbell's 
theorem in Section II-B1 and the aggregated interference model in (8), any two stationary point 
processes with the same intensity have equal mean power of interferences (when the expectation 
of (8) exists). 

Assuming the effective channel power between the desired receiver and its transmitter is h s 




(8) 



g(\\x\\) = {a + b\\x\\>) 



S\-l 



(9) 
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which can be a non-negative constant or a RV, the signal to interference ratio (SIR) is given by 

h s 

SIR = -j (10) 

where / is the interference power given by (8). If there is additive noise with power W which 
can be a non-negative constant or RV as well, then the signal to interference plus noise ratio 
(SINR) is given by [4], [5], [22] 

SINR = WT1- 

In this paper, the LT ordering of two interference distributions will be mainly discussed. The 
Laplace transform of interference plus noise is E[exp (— s(W + 1))] = Cw(s)Ci(s). Clearly 
the LT ordering of two interference distributions is not affected by a common noise power W. 
Therefore, we focus on interference distributions hereafter. 

III. Ordering of Outage Probability and Ergodic Capacity Metric 

In this section we introduce performance metrics involving the stochastic ordering of inter- 
ference. Firstly, we study the SIR-based outage probability. It has been shown in [5] that when 
h s is exponentially distributed, the LT ordering of the interference can be related to the outage 
defined in terms of SINR. We generalize this result to a broader class of distributions for the 
effective channel RV h s . 

Theorem 1. Let I\ and 1 2 denote the interferences from point process $1 and $2 respectively. 
Also, let hs be the effective fading channel between the desired receiver and its transmitter 
and has a CCDF Fh s (x) := 1 — Fh s (x), which is a cm. function. Under these assumptions, if 
h <Lt h, then SIR 2 < st SIRl 

Proof: Let F hs (x) be the cm. CCDF of h s . Then we have 

PiSIR, >x)= E h [P(h s > xh)] = E h [F hs (xh)] (12) 
> E /2 [F hs (xl 2 )] = E h [P(h s > xl 2 )\ = P{SIR 2 > x), (13) 



March 3, 2013 



DRAFT 



8 



where the inequality is due to I x < Lt I 2 and equation (3). Recalling the definition of < st in (1), 



Theorem 1 enables us to conclude the usual stochastic ordering of SIRs whenever LT ordering 
of interferences are established. In the sequel, we will show many examples of how LT ordering 
between two interference distributions can be established which can be used with Theorem 1 to 
establish the "usual stochastic" ordering of outage as in (1). Illustrations of Theorem 1 are shown 
in our numerical results in Section VII. Note that if two interferences have the usual stochastic 
ordering h < st h (a stronger assumption than in Theorem 1), then we can have SIR 2 < s t SIRi 
with any arbitrary distribution for h s and not just those with a cm. CCDF. 

Now consider the expression 



which is the bandwidth-normalized capacity in the weak interference regime [5], [22], [23]. The 
expectation in (14) is with respect to the arbitrary positive random variable hs, as well as the 
interference / given in (8). The noise power W > is fixed, and hs, I are independent. 

We would like to compare two regimes (hs 1 ,h) with capacity C\, and (hs 2 , h) with capacity 
C 2 , where h S2 < Lt hs x , and I-y < Lt I 2 . Since log(l + h s ) is c.m.d. with respect to hs and 
log(l + l/(W+ /)) is cm. with respect to /, by using (3) and (4) successively, we have C\ > C 2 . 
Therefore the LT ordering of interferences and the LT ordering of effective channels together 
lead to ordering of ergodic capacities. Unlike in Theorem 1, in this case, we can compare the 
ergodic capacity regardless of the distribution of fading channel, hs. Section IV-A has examples 
of fading distributions that are LT ordered. 

Having emphasized the impact of LT ordering on outage and capacity metrics, in what follows, 
we will investigate the conditions for the presence of LT ordering of interference distributions. 
The three factors which affect the interference distribution are the fading channel from the 
interfering nodes to the receiver hf\ the path-loss model <?(•), and the point process of interfering 



Theorem 1 is proved. 




(14) 
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nodes $ as shown in the interference model in (8). We derive the conditions on fading channels 
and path-loss models and the underlying point processes for the LT ordering of interferences. 
We also identify the LT ordering of interference distributions in commonly used point processes. 

IV. Comparison of Fading Channels on the Interference Link 

In the previous section we mentioned that the effective channels between the transmitter and 
receiver, h Sl and h s , 2 , can be compared. In this section, we find the conditions on the distribution 

(x) (x) 

of h) on the interference link for the interference in (8) to be LT ordered. Since h) are assumed 
as i.i.d., we will drop the node index x for convenience hereafter, when we refer to its distribution. 

Theorem 2. Let I\ and I2 denote the interferences both with the same path-loss model g{\\x\\), 
from a stationary point process $ as in (8). Also, let and hj 2 be RVs whose distributions 
capture the fading channels between the receiver and interferers under the two scenarios that 
are compared. Under these assumptions, if <u hi 2 , then l\ <Lt h- 

Proof: The proof is shown in Appendix A. ■ 
The interference in a stationary point process depends on the fading channel between the 
receiver and interferers. Intuitively a bigger LT ordering indicates a fading channel that is more 
like AWGN (i.e., "less fading"). Therefore, the interference arising from it indicates interference 
adding up more coherently, giving rise to a "bigger" interference in the LT ordering sense. 

As mentioned in Section II-C, we assumed the interferers are distributed on an infinite region 
due to the stationarity of the point process. However, receivers in many wireless networks may 
experience interference from finite-area regions [24]. Theorem 2 is also valid for the interference 
in finite area. Moreover, when guard zones around the receiver occur due to sophisticated MAC 
protocols [25], [26], Theorem 2 still holds. 



March 3, 2013 



DRAFT 



10 

A. Parametric Fading on the Interference Link 

In this section, we show the interference distributions are monotonic in line of sight (LoS) 
parameter of the fading channels with respect to the LT order for commonly used parametric 
fading distributions such as Nakagami-m and Ricean fading. By Theorem 2 this implies LT 
ordering of interferences, which, by Theorem 1 establishes stochastic ordering of SIRs. 

Consider Nakagami fading model, where the envelope \fh~i is Nakagami and the effective 
channel hi follows the distribution 

f hl (x) = — — x m - 1 exp(-ma;), x > 0. (15) 
T(m) 

Since E[— shi] = (1 + s/m)~ m is a decreasing function of m for each s, it follows that if the 
m parameters of two channel distributions satisfy mi < m 2 then, <u hj 2 where h Il and hi 2 
have normalized Gamma distributions with parameter mi and m 2 respectively. From Theorem 
2, it follows that Ii < Lt I 2 , if mi < m 2 . 

Similarly, the envelope ^fhi is Ricean and the distribution of effective channel hi is given by 

f hl (x) = (K+1) exp [—{K + l)x - K\ I (2^K(K + l)x), x > 0, (16) 

where K is the LoS parameter of Ricean fading channel and I (t) := J2m =0 (t/2) 2m /(m\T(m + 
1)) is the modified Bessel function of the first kind of order zero. The Laplace transform of (16) 
decreases with K for all s > 0. Thus, similar to the Nakagami case, if K 1 < K 2 are the Ricean 
parameters of two channels, then h Il < Lt h l2 which by Theorem 2 imply that l x < Lt I 2 . 

In addition to the Nakagami-m and Ricean fading cases, it can be shown through a procedure 
similar to the discussion that the interference distribution corresponding to Nakagami-g (Hoyt) 
fading [27] also satisfies LT ordering with respect to the shape parameter. 

B. Interference in Combined Multipath Fading and Shadowing 

The effect of shadow fading on the interference power distribution can be modeled as a product 
of a shadowing random variable with a multipath fading random variable. Let h^ <u hj 2 be the 
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two effective multipath fading distributions, and Xi < Lt X 2 be the two shadowing distributions. 
Then, from [10, Theorem 5.A.7 (d)], it follows that the composite RV satisfies hi 1 X 1 < Lt h l2 X 2 , 
since l(x,y) = xy has a cm. derivative in each variable. We conclude that if < Lt h l2 and 
X\ <Lt X 2 , then Ji < Lt J 2 from Theorem 2. Such conclusions are especially useful even when 
the composite distribution of h^Xi or hj 2 X 2 cannot be written in closed-form. 

V. Comparison of Path-loss Models 

Here, we show the ordering of mean power of interferences in a stationary point process with 
different path-loss models. Generally, the mean power of interference is an important factor to 
determine a system performance. For example, by considering the mean power of interferences as 
noise power, the average error rate can be approximated. Thus, we are interested in comparison 
the mean power of interferences. In what follows, we show the condition on the mean power of 
interferences with non-singular path-loss models with different path-loss exponents under which 
this ordering holds: 

Theorem 3. Let I\ and I 2 denote the interferences both with identical fading distribution hi in a 
stationary point process $. Also, let gi(\\x\\) and fi^GI^H) be the non-singular path-loss models 
with a = 1 and b = 1 in (9) and with different path-loss exponents 5\ and 5 2 , respectively. If 
£i < 5 2 , then E[J 2 ] < E[Ji]. 

Proof: The proof is given in Appendix B. ■ 
Theorem 3 compares the mean power of interference for two different path-loss models and 

identical fading. It is clear from the proof of Theorem 3 that if the fading is different with 

E[/i/J > E[/i/ 2 ], the conclusion would still hold. 

Under a stationary Poisson point process, we can establish the stronger stochastic ordering, 

LT ordering of interferences, because if Ii > Lt I 2 , then E[7i] > E[7 2 ]- In what follows, we find 

the conditions on the interference distributions with non-singular path-loss models with different 
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path-loss exponents under which this stochastic ordering holds: 

Theorem 4. Let I\ and I 2 denote the interferences both with identical fading distribution hi 
in a stationary Poisson point process $ppp. Also let ^i(||x||) and g 2 i\\x\\) be the non-singular 
path-loss models with a = 1 and b = 1 in (9) and with different path-loss exponents 5\ and 5 2 , 
respectively. If 5± < 5 2 , then I 2 <u h- 

Proof: The proof is given in Appendix C. ■ 
For the non-singular path-loss model gr( ||x||) with a = 1 and b = 1 is assumed, however, the 
mean power of interference with is greater than that with <?2(||£||) from Theorem 3. 

Thus, for more fair comparison of LT ordering of the interferences, we set the parameter b in 
(9) for g2(||^||) to have a equal mean power of interference as that of using Campbell's 

Theorem as follows: 



a = 1 and 6=1 and g^OI^H) with a = 1 and b is set the value in (17) and with different 
path-loss exponents 5i and 5 2 in (9) in a stationary point process $ have same mean power 
of interference, even if 5i < 5 2 . Then, with (17), we can establish additional LT ordering of 
interferences as follows: 

Corollary 5. Let l\ and l 2 denote the interferences both with identical fading distribution hi in 
a stationary Poisson point process $ppp. Also let <7i(||a;||) be the non-singular path-loss model in 
(9) with 81, a — 1 and 6=1 and g^GI^H) be the non-singular path-loss models with 5 2 , a = 1 
and b is set the value in (17). If Si < 5 2 , then I 2 <u I±. 

Proof: The proof is the same as the proof for Corollary 4 with setting the value of b for 
g 2 (\\x\\) is (17) instead of b = 1. ■ 
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From Corollary 5, it is seen that the interference distributions with non-singular path-loss 
models in stationary Poisson point process are monotonic in the path-loss exponent with respect 
to the LT order. Indeed, the parameter b in (17) was set to ensure the mean interference powers 
are equal. Hence, the better performance in SIR-based outage probability or ergodic capacity 
with increased path-loss exponent 5 is not due to an improvement of average interference power. 

Theorem 3, 4 and Corollary 5 show that one path-loss model need not dominate the other 
pointwise to get an ordering in interference in the case of the non-singular path-loss model. 
Indeed, if g\(r) > g 2 (r),\fr e M + (for example, a comparison of singular and non-singular path- 
loss model), it is obvious that E[Ji] > E[J 2 ] since the aggregated interference power with gi(r) 
is always greater than that with g2(r). For a stationary Poisson point process, we can also easily 
observe when gi(r) > #2 (?"■)> ^ — u h from (27). However, from Theorem 3, 4 and Corollary 
5, we can establish the stochastic orderings of interferences even in the case that gi(r) < g 2 (r) 
in some range of r and gi(r) > g2(r) in another range of r as shown in Fig. 2. The simulation 
result to verify these theorems will be shown in Section VII. 

For the singular path-loss model (a = and b > 0) in (9), a finite mean of interference 
power does not exist in a stationary point process since (23) does not converge. Thus the means 
of interferences cannot be compared. In case of a stationary Poisson point process, we have a 
closed-form expression for Laplace transform of interference with the singular path-loss model 
in a stationary Poisson point process as follows [22]: 

d(s) = exp(-Ac d E[/#]r(l - a)s a ), (18) 

where a = d/8 and a is the volume of the rf-dimensional unit ball. E[hf] is a fractional moment 
of hi with < a < 1. Unlike the non-singular path-loss model, however, the LT ordering does 
not hold between two different interferences corresponding to path-loss exponents 8 1 < S 2 in 
case of the singular path-loss model since (18) is not a decreasing function of 5. 
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VI. Comparison of Different Point Processes 

In [14] directionally convex ordering (dcx) is used to order point processes. In this section, we 
define a new stochastic ordering between point processes and state some results involving the 
ordering of point processes. We first define a new stochastic ordering of point processes based 
on the well-known Laplace functional: 

Definition 1. Let $x and $ 2 be two stationary point processes such that 



L*>Au) = E 





> E 











= W«) (19) 



where u(-) runs over the set U of all non-negative functions on M d . Then $i is said to be smaller 
than $ 2 in the Laplace functional (LF) order (denoted by $i <Lf $2)- 

It can be shown that LF ordering follows from dcx ordering, which makes LF ordering easier 
to verify and easier to relate to interference metrics. Note that the LT ordering in (2) is for RVs, 
whereas the LF ordering in (19) is for point processes. They can be connected in the following 
way: 

$1 <Lf $2^$^ u(x) < u u ( x )> VueU (20) 

Hence, it is possible to think of LF ordering of point processes as the LT ordering of their 
interferences in the absence of fading, for all non-negative path-loss functions. But as we will 
see, LT ordering of their interferences in the presence of fading can be proved when two point 
processes are LF ordered. We next prove a generalization of Theorem 2 where the two point 
processes are different and LF ordered. 

Theorem 5. Let $1 and $ 2 be two stationary point processes and and h l2 be RVs whose 
distributions capture the fading channels between the receiver and interferers under the two 
scenarios that are compared. Also let l\ and I 2 denote the interferences with same path-loss 
model g(\\x\\) in $! and <3> 2 respectively. If $1 <Lf $2 and <Lt hj 2 , then Ii < Lt I 2 . 
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Proof: Theorem 5 follows from Theorem 2 and equation (20) with u(x) = hig(\\x\\). ■ 
Note that in Theorem 5, hi does not only capture a fading distribution, but also can capture 
a random thinning property of point process. 

A. Laplace Functional Ordering between Specific Point Processes 

Here, we compare the LF order between different point processes which are commonly used. 
These are Poisson point process and Neyman-Scott process with Poisson distributed number of 
points in each cluster which is one of example of Poisson cluster process. The Neyman-Scott 
process results from homogeneous independent clustering applied to a stationary Poisson point 
process. The details for these point processes can be found in [3], [16]. 

Theorem 6. If $ppp and $pcp denote Poisson point process and Neyman-Scott process with 
Poisson distributed number of daughter points respectively and both point processes have same 
intensity \, then $ PC p <Lf $ppp- 

Proof: The proof is given in Appendix D. ■ 
We consider another point process, the mixed Poisson process which is a simple instance 
of a Cox process. It can be thought of as a stationary Poisson point process with randomized 
intensity parameter X which has the averaged intensity measure Ex[X] = A. Every sample of 
such a process looks like a sample of some stationary Poisson point process. We compare the 
LF ordering of this point process with that of stationary Poisson point process as follows: 

Theorem 7. If $ PPP and $mpp denote Poisson point process and mixed Poisson process respec- 
tively and both point processes have same average intensity \, then <J>mpp <Lf ^ppp- 

Proof: The proof is given in Appendix E. ■ 
The aggregated interference from Poisson cluster process or mixed Poisson process is always 
less than that from Poisson point process in LT order by Theorem 5. This means that the 
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orderings of SIR-based outage probabilities or ergodic capacities are established depending on a 
point process by the performance metrics in Section III. It is noted that from Campbell's theorem 
in Section II-B1, any two stationary point processes with same intensity have equal mean power 
of interferences (when the expectation in (8) exists). Hence, the better performance in SIR-based 
outage probability or ergodic capacity in specific point process is not due to an improvement of 
average interference power. 

The above point processes are less than Poisson point process in LF order. In what follows, we 
show the point process which has larger LF ordering than Poisson point process even though it is 
non-stationary and non-isotropic. In binomial point processes, there are a total of N transmitting 
nodes uniformly distributed in a rf-dimensional ball of radius r centered at the origin, denoted 
as B (r). The density of the process is given by A = N/(cdr d ) where Cd is the volume of the 
d-dimensional unit ball [28]. In this case, we can compare the LF ordering of point processes 
in bounded area as follows: 

Theorem 8. Let $ppp(r) be a Poisson point process over B (r). //" $ppp(r) and $bpp denote Pois- 
son point process and binomial point process respectively infinite area and both point processes 
have same intensity A, then $ppp(r) <Lf $bpp whenever < A f B() ^ [1 — exp(—u(x))] dx < N 
holds. 

Proof: The proof is given in Appendix F. ■ 
In a finite area B (r), the aggregated interference at the origin from a binomial point process 
is always larger than that from Poisson point process in LT order from Theorem 5. 

B. Laplace Transform Ordering of Interferences in Heterogeneous Networks 

As unlicensed band utilization increases, the unlicensed wireless network may experience 
adverse interference from collocated wireless devices that are transmitting in the same unlicensed 
band. Such a heterogeneous network scenario can be modeled as a superposition of mutually 
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independent point processes [29]. Let $1 = (J^ <E» 1; j and $ 2 = $ 2 ,j for i = 1,...,M 
be the heterogeneous networks which are modeled as superpositions of mutually independent 
point processes. Since the Laplace functional of superposition of mutually independent point 
processes is L$(u) = YlfLi L ^i( u ) [30], if <Lf $2,* for i = 1,...,M, then $i < Lf $ 2 . 
Therefore, 7i <Lt ^2 from Theorem 5. 

VII. Numerical Results 

In this section, we verify our theoretical results through Monte Carlo simulations. Since the 
LT ordering between two interference scenarios cannot be verified directly from its probability 
distributions such as PDF and CDF, we will verify the LT ordering between interferences by the 
ordering of SIR-based outage probabilities or ergodic capacities as mentioned in Section III. 

A. Comparison of Fading Channels on the Interference Link 

In many practical scenarios, different links in wireless networks can experience asymmetric 
fading conditions. If the interferer's channel is Nakagami-m fading, while the desired link is 
Rayleigh fading which has an effective channel CCDF which is cm., we can compare SIR-based 
outage probabilities using using Theorem 1. 

In Fig. 3 the CDFs of interference power and SIR from Poisson cluster process with different 
Nakagami-m fading parameters and with the non-singular path-loss model (a = 1, b = 1 and 
5 = 4) in (9) are shown. We consider two different LoS parameters: m x = 1 and m 2 = 2. 
The choice of these parameters ensures l\ < Lt I 2 from Theorem 2. Consequently, we observe 
SIRi >st SIR 2 in the bottom of Fig. 3 which agrees with Theorem 1 even though there is a 
crossover point between interference power distributions in the top of Fig. 3. 

The CDFs of SIR from Poisson point process with different Nakagami-m fading parameters 
and with the singular path-loss model (a = 0, b = 1 and 5 = 4) in (9) are shown in Fig. 4. 
Similarly the LoS parameters mi < m 2 lead to I\ <Lt h from Theorem 2. Clearly, it is observed 
SIRi >st SIR 2 in Fig. 4 which agrees with Theorem 1. 
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Table I shows the ergodic capacity performances when the desired link has Ricean fading 
channel with K s = 5 and interfering channels follow Ricean distributions with K Il = and 
Kj 2 = 1 in the two scenarios compared on a Poisson cluster process. The ergodic capacity with 
Ii is always better than that with I 2 as expected since the interference distributions are monotonic 
in LoS parameter of Ricean fading channel with respect to the LT ordering in Section IV. 

The ergodic capacities in Poisson point process are shown in Table II. In this case, all 
conditions are same except for the type of point process. Therefore, the ergodic capacity with 
7i is always better than that with I 2 when h < Lt I 2 . 

B. Comparison of Path-loss Models 

We show in Fig. 5 the CDFs of the interference power and CDFs of SIR from a Poisson point 
process with the non-singular path-loss models, gi(||a;||) and p 2 (||^||) which are given in Fig. 2 
and discussed in Corollary 5. It is noted the non-singular path-loss models with two different 
path-loss exponents: S± — 4 and 5 2 = 8 ensures h >u h from Corollary 5. We consider 
additional non-singular path-loss, g 3 (||:r||) whose parameters are a = 1, b = 1, and 5 = 8 in Fig. 
2. Since S^dl^H) < fi^GI^H) for ||x|| > as shown in Fig. 2, it is obvious I 2 > Lt 7 3 . Thus, 
we observe SIRi < st SIR 2 < st SIR 3 in the bottom of Fig 5. By Theorem 1, when the fading 
channel between the desired receiver and its transmitter is Rayleigh distributed and the effective 
fading channel is exponentially distributed, we can observe the usual stochastic ordering of SIR 
distributions if the interferences are LT ordered as shown in Fig. 5. 

C. Comparison of Different Point Processes 

In the following, we compare the CDFs of interference power and CDFs of SIR from Poisson 
point process and Poisson cluster process. Since the interferences from Poisson point process 
and Poisson cluster process with the same intensity A, the same non-singular path-loss model 
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g(||x||), and identical fading distribution hi yield / PCP < Lt / P p P , which implies SIR PCP >st SIR PPP . 
This is observed in the bottom of Fig. 6 as predicted from our theoretical result in Theorem 6. 

In addition to SIR-based outage performance, it is observed that the ergodic capacity in Poisson 
cluster processes is always greater than that in Poisson point processes by comparing same rows 
in Table I and II. 



In this paper, we used stochastic orders to compare performance in wireless networks. We 
showed that when interference is LT ordered, it is possible to order the SIR in the usual stochastic 
ordering sense when the effective channel has a cm. CCDF. Similar results hold when the metric 
is the bandwidth-normalized capacity. This lead to the study of the conditions for LT ordering 
of interference. Three factors affecting interference are the fading channel from the interfering 
nodes to the receiver, the path-loss model and the distribution of the interfering node location. We 
derived conditions on these factors so that LT ordering between interferences holds. In addition, 
we defined Laplace functional ordering of point processes and derived its inherent stochastic 
ordering of interferences when the fading channel and the path-loss model are assumed to be 
same for both point processes. The power of this approach is that such comparisons can be made 
even in cases where a closed form expression for the interference is not analytically tractable. 
We verified our results through Monte Carlo simulations. 



The Laplace transform of aggregated interference / := J2 x e^^i 9(\\ x \\) * s a Laplace func- 
tional (6) evaluated at u(x) = shig{\\x\\), s > 0, hi > and gf(||x||) > where hi is the effective 
fading channel between the receiver and interferers and is a path-loss model. £/(s) can 

be expressed as follows: 



VIII. Summary 



Appendix A: Proof of Theorem 2 



Ci(s) = E [e- sI ] = E e 



shi f Rd g(\\x\\)$(dx) 



(21) 
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where the expectation is to be taken over both $ and hj. Let Z = J Rd g(\\x\\)$(dx) in (21). 
From [10, Theorem 5.A.7 (b)], if E [exp (sh)\Z = z] > E [exp (sI 2 )\Z = z] for all z in the 
support of Z, then I x < Lt J 2 - Therefore, it is sufficient to show the following equation regardless 
of a point process $ in order to satisfy the LT ordering between interferences, 

E h/i [e- szh 'i \Z = z}> E hh [e- szh <*\Z = z] . (22) 

But (22) follows from the assumption <u h l2 . Thus, we conclude that if hj x < Lt h l2 , then 

h <Lt ^2- 

Appendix B: Proof of Theorem 3 

In order to prove E[7i] > E[7 2 ], using Campbell's theorem in Section II-B1 with u(x) = 
hig(\\x\\) where hi is the (power) fading coefficient and ^(||x||) is the non-singular path-loss 
model in (9), we need to show the following: 

\[ hi 9l (\\x\\)dx > X [ hig 2 (\\x\\)dx. (23) 

JR d JR d 

Since hi is independent from the point process and can be expressed as g(r),r = \\x\\ 

under polar coordinates, after expectation with respect to hi and change to polar coordinates, 
the following condition needs to be satisfied to prove Theorem 3: 

E[h T }Xc d d / gi (r)r d - l dr > E[hi] Xc d d \ g 2 (r)r d ~ l dr. (24) 
Jo Jo 

where c d is the volume of the d-dimensional unit ball and hi is the (power) fading coefficient 
between the receiver and interferers. Using a change of variables, we get 

/ u (^r 1 ^)) 7T- {9i\n)) du> u {g 2 \u)) d - 1 — {g 2 \u)) du. (25) 

J gi (oo) ° U Jg 2 (oc) OU 

where d is the dimension of point process, is the inverse function of g(-) and I[u e S] — 1, 

if u e S, and otherwise, is the indicator function. Substituting the non-singular path-loss models 
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gi(r) and g 2 (r) with a = 1 and 6 = 1 in (9) into (25), we get 




S 2 u 2 J 5i 5 2 

(26) 

since " r 5 f ys> is a decreasing function with 5 for a fixed d and 5i < S 2 . The proof for 
Theorem 3 is complete. 

Appendix C: Proof of Theorem 4 

The Laplace transform of interference power in a stationary Poisson point process with path- 
loss model gj(r),j = 1,2 in (9) and intensity A can be expressed as follows: 

= exp {"Mid jf t 1 - ex P (shigj(r))] r^Mrj (27) 

where q is the volume of the <i-dimensional unit ball and hj is the (power) fading coefficient 
between the receiver and interferers. From [10, Theorem 5.A.7 (b)], it is sufficient to show 
£h(s) > £/ 2 ( s ) m (27) regardless of a distribution of h I in order to satisfy the LT ordering 
between interferences. To do so, the following condition needs to be satisfied after change of 
variables: 

/ [1 - exp (shju)} (g^H) (9T\u)) du 

J 9l (oo) OU 

1-92(0) a 

> / [l-e W (-sh I u)](g; 1 (u)) d ' — (g^iu^du. (28) 

Jg 2 (oo) CU 

Substituting the same non-singular path-loss models into (28), it follows 



Jo 



1 " -i d i 

i i-i ^ 



1 - exp (shju)) I vu _L | dM (29 ) 



> (1 - exp (sh:)) I u I i*L_J. U ^ | d M (30) 

(1 - exp (-**,)) | r ( 1 ~^) r (^) _ r ( 1 ~g) r U/ | > Q _ (31) 
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(30) follows from (1 — exp(— c))u < 1 — exp(— cu) for c > and < u < 1 and (31) follows 
from 1 — exp(— shi) > for s/i/ > and A > from (26). Theorem 4 is proved. 

Appendix D: Proof of Theorem 6 

Let u(-) G U be all non-negative functions on M. d . The Laplace functional of Neyman-Scott 
process with Poisson distributed number of daughter points and with the distribution /(•) for 
locations of daughter points can be expressed as follows [3], [5]: 

L^ CP (u) = expj-Ap^ 1 - exp^-c^l - J ^exp(-u(x + y))f(y)dyj^ dxj 



> exp\-\ p c 1-/ exp(-u(x + y))f(y)dy 
= exp < —A / [1 — exp(— u(x))] dx > 



da; 



(32) 



(33) 



— L$ m (u) 

where the inequality in (32) follows from the fact that 1 — exp(— ax) < ax, a > and (33) 
follows from change of variables, interchanging integrals and using J f(y)dy = 1 [16]. 

Appendix E: Proof of Theorem 7 

Let u(-) E U be all non-negative functions on R d . The Laplace functional of the mixed Poisson 
process process with a random intensity measure X which has the averaged intensity measure 
Ex[^] = A can be expressed as follows [31]: 



£*MP P («) = E X 



exp < —X / [1 — exp(— u(x))} dx > 

I JR d ' J 

> exp \ -E X [X] / [1 - exp(-u(a;))] dx 

I jR d 

= exp < —A / [1 — exp(—u(x))] dx > 

I JR d J 



(34) 
(35) 
(36) 



where the inequality in (35) follows from Jensen's inequality since the term inside the brackets 
in (34) is a convex function of X. 



March 3, 2013 



DRAFT 



23 



Appendix F: Proof of Theorem 8 

Let u(-) G U be all non-negative functions on R d . The Laplace functional of the binomial 
point process consisting N points with a density A can be expressed as follows [28]: 

L* BPP (u) = (l- Ajf [l-exp(-u(x))]dr) (37) 
< exp < —A / [1 — exp(—u(x))] dx > 

I JB (r) J 
= L $ppp(r) (M) 



(38) 



where the inequality in (38) is due to (1 — c/n) n < e c for < c < n. Thus, Theorem 8 is 
followed whenever < A j Bo ^ [1 — exp(—u(x))] dx < N holds. 
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TABLE I 

Ergodic capacities (bits/s/Hz) over Ricean fading channel with Ks = 5 in Poisson cluster process 



SINR (dB) 


-4 


-2 
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4 


6 


8 


10 


h (K h = 0) 


0.9433 


1.2282 


1.5709 


1.9820 


2.4538 


3.0063 


3.6482 


4.5313 


h (K l2 = 1) 


0.9426 


1.2277 


1.5707 


1.9803 


2.4526 


3.0020 


3.6449 


4.5278 



TABLE II 

Ergodic capacities (bits/s/Hz) over Ricean fading channel with Ks = 5 in Poisson point process 



SINR (dB) 
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h (K h = 0) 


0.9346 


1.2162 


1.5485 


1.9501 


2.3847 


2.9295 


3.5491 


4.3477 


h (K l2 = 1) 


0.9342 


1.2152 


1.5468 


1.9460 


2.3816 


2.9231 


3.5407 


4.3349 
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Fig. 1. Illustration of a wireless network. The black dots represent interfering nodes which form a point process $ and the 
dotted lines represent their interfering signals. The white dot and the triangle at the origin are the desired transmit/receiver pair 
which are not part of the point process. 



Non-singular path loss models 



g 1 (||x||)w/8=4,a=1,b=1 




Fig. 2. Non-singular path-loss models with different path-loss exponents, Si and 82 
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CDF of interference power 




Fig. 3. CDFs of interference and SIR for Poisson cluster process with different fading parameters and with A = 0.01 
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Fig. 5. CDFs of interference and SIR for Poisson point process with non-singular path-loss models with different path-loss 
exponents and with A = 0.01 



CDF of interference power 




CDF of SIR 

0.3 1 1 1 




Fig. 6. CDFs of interference and SIR for Poisson point process and Poisson cluster process with A = 0.01 
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